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(54) Method and apparatus for determining the concentration of a component 



(57) A method of determining at least the concen- 
tration of a component from the intensity of electromag- 
netic waves with at least two selected wavelengths 
which are reflected by human tissue a transmitted 
through human tissue comprises firstly the step of con- 
verting the intensities of the received electromagnetic 
signals into at least one first and one second time- 
dependent electric signal (R it IRJ. Then a time-discrete 
transformation of the first and of the second electric sig- 
nal into the frequency domain is performed to determine 
first and second spectral values (Sr*. Sirk) of the first 
and of the second signal. Complex combinatorial values 
(A*, ratio*) are formed from said first and second spec- 
tral values and physiologically relevant combinatorial 
values are selected by evaluating the complex combina- 
torial values according to given criteria for the physio- 
logical relevance thereof. Finally, the concentration of 
the component is calculated by using the selected com- 
binatorial values or by using the frequencywise-associ- 
ated spectral values. 
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dependent electric signal; 

timeKliscretely transforming the first and of the second electric signal into the frequency domain for determining 
first and second spectral values of the first and of the second signal; 

forming complex combinatorial values from said first and second spectral values; 

selecting physiologically relevant combinatorial values by evaluating the complex combinatorial values according 
to given criteria for the physiological relevance thereof; and 

calculating the concentration of the component making use of the selected combinatorial values or making use of 
frequencywise-associated spectral values. 

According to a preferred embodiment of the present invention, for determining the correlation ol a cornponent 

** r^SLZ SSJuHV. .uW* time window, in acoordane. wilh a prefe-nrf .MM 

ZZZZiZZZ n-St coo*ir*lonal values H oamed out on th. M<^>«>^» 

^XZ^mZw** £ ""W con*H«orial valuos. In accorfanoe wHn preferred enMm*** the ixewrt 
;^r^.*^r,l and .%Mta. ind« are dewmm*. on the ba* of the seteeted e»*™w»l »W <* 

r^taadtnTe^ TofaSeXctei can be subjected Bfite,ing^™raoln B .™egi^cnw« U sad cana**- 
lowing step6: 

converting the intensities of the received electromagnetic signals into at least one first and one second time- 
dependent electric signal; 

time-discretely transforming the first and of the second electric signal into the frequency domain for determining 
first and second spectral values of the first and of the second signal: 

selecting physiologically relevant first and second spectra, values * evaluating said first and second spectra, val- 
ues according to given criteria for the physiological relevance thereof, and 
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calculating the concentration of the component making use of the selected first and second spectral values. 

The present invention additionally provides apparatus for carrying out the methods described hereinbefore. 
Further developments of the present invention are disclosed in the dependent claims. 
5 In the following, preferred embodiments of the present invention are explained in detail making reference to the 

drawings enclosed, in which: 

Fig. 1 shows a schematic block diagram representing an oximeter: 
10 Rg. 2 shows a representation of preprocessed raw signals I and IR in the time domain; 

Fig. 3 shows a representation of the Fourier transformed values of the signals shown in Fig. 2; 
Fig. 4 shows a complex representation of the signals shown in Fig. 3; 

15 

Fig. 5 shows a representation of the amount spectrum of complex combinatorial values determined on the basis of 
the signals shown in Fig. 3; and 

Fig. 6 shows an overview of a preferred embodiment of the method according to the present invention in the form 
20 of a flow chart 

Rg. 1 shows the basic structure of an oximeter which is adapted to be used for carrying out the method according 
to the present invention. The oximeter includes a sensor device 1 0 for transmitting and receiving electromagnetic waves 
with at least two selected wavelengths. In accordance with a preferred embodiment of the present invention, one wave- 

25 length in the visible red spectrum and another wavelength in the infrared spectrum are used for this purpose. The signal 
produced by the photoelectric receiver in the sensor device 10 is fed through a line 12 into a current-to-voltage con- 
verter unit 14 and. subsequently, into an analog-to-digrtal converter 16 producing a digital representation of the intensity 
measured. This digital value is then fed into a microprocessor 18 working under the control of a program memory 20. 
The program memory 20 contains the whole code which is necessary for the processor for determining the oxygen sat- 

so uration and, if desired, also a perfusion index or a pulse rate, as will be explained in detail hereinbelow. The program 
code in said memory 20 instructs the processor 1 8 to carry out all the steps which are necessary for the method accord- 
ing to the present invention. 

It is obvious that, instead of the special arrangement described hereinbefore, ft is possible to use arbitrary process- 
ing units which are known in the field of technology and which are capable of carrying out the methods according to the 
35 present invention. Such processing units can, for example, comprise display and input devices and additional memory 
and accessory arrays in the manner known. 

According to preferred embodiments of the present invention, the electric raw signals, which have been obtained 
on the basis of the intensity of electromagnetic waves at at least two selected wavelengths, typically red and infrared, 
are first subjected to preprocessing so as to remove time-dependent drifts from said first and second signals. In Rg. 2, 
40 raw curves R (red) and IR (infrared), which have been subjected to such preprocessing, are shown. 

In the following, a brief explanation of this preprocessing will be given. The preprocessing aims at separating, if 
possible, all the time-dependent superpositions from the raw signals. For this purpose, a continuous average value is 
formed on the basis of the raw waves Rj and IRj for determining a so-called baseline according to the following equa- 
tions: 

45 

hT/2 

M ro- 1/r 2 R i (3a) 

SO UT& 

M |Ri «1/T £ 'R| (3b) 

55 wherein i = the continuous sample; 

Mr,. M (Rj = the continuous average values of Rj and IRj symmetrically around said sample i; 
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j = period of averaging (e.g. one second). 

The baseline is now eliminated from the raw signals by means of the following equations: 

6 R , = ln(R,/M Rl ) W 

s IRl »ln(IR,/M IRI ) (4b) 

s Ri and s IRi now represent the continues samples of the original raw signals for red and infrared in the time domain. 
Examples of such time-dependent samples s Ri . curve 22. and s, Ri . curve 24. are shown in the form of curves in Fig. 
2 The ordinate represents the time in the diagram according to Fig. 2, whereas the abscissa represents normalized 
amolitudes of s Rj and s IRi . In view of the fact that the raw curves were sampled at 1 25 Hz in the above example, but are 
strongly band-limited to 10 Hz. redundancy was additionally eliminated in the course of the preprocessing and sampled 

down by the factor of 4 to 31 .25 Hz. . 

The above-described preprocessing of the raw signal provides a plurality of advantages. On the one hand, the sig- 
nals are normalized, whereby the range of values becomes smaller and this is advantageous with regard to integer 
processing Furthermore, the signal ratio of each sample pair R to IR already provides a ratio that can be utilized 
according to equation (2). Furthermore, linear drifts of the signals are completely eliminated by the symmetric averag- 

According to the present invention, the signals present in the time domain are now subjected to a time-discrete 
transformation into the frequency domain. According to the preferred embodiment, the preprocessed samples e ra and 
s are subjected to a fast Fourier transformation (FFT) making use of suitable windowing. The selection of the window- 
shape is uncritical in this connection. In the embodiment described, the known cos-shaped Hanning window has been 
used which combines advantageous properties with regard to a minimum peak enlargement and small secondary 
2s Deaks What is more important is the window length. In this connection, a compromise must be made between a win- 
dow length that is sufficiently long for realizing a good frequency resolution, i.e. many points, and a window length that 
is sufficiently short for rapid changes of the signal frequencies, of the pulse rate and of intermittent disturbances. In the 
oreferred embodiment a window of 8 seconds, n = 256 samples, was used. Making use of such a window, a sufficient 
frequency resolution of 1/(32ms • n) = 0.12 Hz is realized. Hence. 82 values are obtained in the frequency band of 

so i nter ^ t e fr F ^°eI 0 trlnsformed values Sr*. curve 26. and S IRK . curve 28. of the signals of Fig. 2 are shown in Rg. 3. said 
Fourier transformed values being determined in the way described hereinbefore. The ordinate of the diagram of Fig 3 
represents the frequency, whereas the abscissa again represents normalized amplitudes of Sr* and Sir*. Due to the 
baseline preprocessing described hereinbefore, the DC component is practically equal to zero with the exception of 

35 '"''when'?? def ined as amplitude spectrum (absolut value of the Fourier transformed values) of the time function s 
(equations 4a and 4b). a ratio can be formed on the basis of the ratio of the coefficients for each pair of frequency 
points: 

Ratio h-SRu/S,™, © 

wherein k stands for the respective frequency base points. 

For the peaks of the spectrum, which originate from a Wood pulse, a ratio should be obtained in this way which 
leads to the true Sp02 value (i.e. oxygen saturation value). For spectral components lying outside of the regions around 
the fundamental and harmonic waves of the Wood pulse, ratio results are obtained which must be regarded as disturb- 
ers and which must be eliminated. This elimination will be explained in detail hereinbelow. 

In a continuous mode of operation, the FFT window is advanced by a predetermined period of time. e.g. 1 second, 
in each case. Hence, a new spectrum (pair) is determined for each predetermined period of time. Although this results 
in a 7/8 window overlap and. consequently, in great redundancy, it turned out that the spectral changes occurring in the 
course of 1/8 displacement already justify a new calculation cycle. The method according to the presenUnvenliQiLis 
therefore preferably run through cyclically, the duration of one cyde being e.g. 1 second. 

When in the case of an alternative embodiment of the method according to the present invention, the above- 
described preprocessing of the time-dependent signals is dispensed with, the baseline or the reference value can be 
determined in the frequency domain in some other manner. e.g. as the DC component of the signal m the FFT window 
or as output signal of a normal causal lowpass f Oter having a sufficiently low limiting frequency (e.g. in the region of 0.5 
Hz) In this case the fast Fourier transformation is carried out directly on the raw signals obtained. 

In the preferred embodiment of the method according to the present invention, complex combinatorial values are 
now formed on the basis of the Fourier transformed values. In order to make this more easily comprehensible, it can be 
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considered as a parametric x-y representation of the Fourier transformed values. When the infrared spectrum is plotted 
in the x -direction and the red spectrum in the y-direction, a representation is obtained which has needlelike tips. Such 
a representation for the spectra shown in Fig. 3 is shown in Fig. 4. These "needles" correspond to the peaks of the spec- 
tra. For undisturbed signals very slim needles are obtained, the respective needles for the fundamental and harmonic 

5 waves lying one on top of the other. The direction of the needles corresponds to the saturation. This means, in concrete 
terms, that the mean gradient of a needle is the ratio searched lor. 

Various disturbance frequency components lie outside of the needle generated by a blood pulse. Background noise 
and minor disturbances cause an accumulation of points around the origin in the x-y representation, which could be 
referred to as "cloud". Furthermore, said background noise and said minor disturbances cause an offset of the needles. 

10 Also correlated disturbances of the red spectrum and of the infrared spectrum, such as motion arti f ac t s, which have 
similar spectral components can be discerned as extra needles and are normally not in the regression with the useful- 
signal needles of the blood pulse spectrum. The broken line 40 in Fig. 4 indicates an SpQ2 Emit of 0 percent, whereas 
the broken One 60 indicates an Sp02 limit of 100 percent 

For algorrthmically identifying the "needles", it is first of all necessary to determine the length of said needles. For 

75 this purpose, a so-called distance spectrum is calculated on the basis of the amplitudes § and S ^ of the red spec- 
trum and of the infrared spectrum, the values of said distance spectrum being determined according to the following 
equation at each frequency k: 

A k - aJs&k + sfrat (6) 

20 

These values A* therefore represent the amount of the complex combinatorial values determined on the basts of 
the red and infrared spectra. The phase of the complex combinatorial values is given by equation (5). The distance 
spectrum, i.e. the amount spectrum of the complex combinatorial values, is shown in the form of a line chart in Fig. 5. 
Sections 70, 72, 74, 76, which are identified as peak area (= needle) in the future identification, are shown in said Fig. 5. 

25 For identifying "needles", it is now necessary to search the distance spectrum algorrthmically for peaks, i.e. for 
maxima and the associated foot points. Such an algorithm must then identify from the distance spectrum the peaks 
which belong to the undisturbed Pleth wave and which supply the correct values, consequently. Optionally, the algo- 
rithm should additionally permit a certain smoothing of the outputted curves by an elimination of outliers, e.g. by means 
of a median filter and, possibly, by additional temporal averaging. Furthermore, the algorithm should preferably provide 

30 a special indication, if the determination of a saturation value is not possible on the basis of the data available. 

The distance spectrum must now be examined so as to determine peak areas associated with the so-called "nee- 
dles", cf. Fig. 4. In the following, a preferred embodiment for identifying peak areas in the distance spectrum wfll be 
explained; for the sake of simplicity, these areas will only be referred to as peaks hereinbelow. It is, however, pointed 
out that the embodiment described is only one preferred embodiment for identifying peaks; alternative methods, which 

55 are suitable for identifying peaks and which deviate from the method described, can be used as well. 
In the preferred embodiment for indentifying peak areas, the following steps are carried out 

1 . search for highest peak that fulfills the maximum criterion, cf. hereinbelow; 



40 2. search foot point on the right that fulfills the foot point criteria, cf. hereinbelow; 

3. search foot point on the left; 

4. eliminate the found peak between the foot points from the spectrum; 

45 

5. repeat steps 1 to 5 until a maximum number of 10 peaks has been found or the residual spectral lines are below 
a minimum threshold, e.g. 3%; and 

6. accept all peaks that fulfil the peak criteria, cf. hereinbelow. 

50 

Maximum Criterion 

The maximum does not lie at the edge of the spectrum, i.e. for the contact edges, i.e. the foot points, of a cut-out 
peak 10 Hz are the upper limiting frequency and 0.5 Hz are the lower limiting frequency. 

55 

Foot Point Criteria 

A foot point must futf il one of the folowing criteria: 
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- it must abut on a peak which has already been eliminated, but at least at a distance of one line from said peak; 
new foot point = old peak edge. Hence, a common foot point exists; 

_ the amplitude of two successive lines lies below 25% of the amplitude of the peak ==> foot point = first line; 

- the amplitude of two successive lines lies below 50% of the amplitude of the peak and said lines descend monot- 
onically and gently, i.e. the change in height normalized to the peak height is < 5% ==> foot point = line in front of 
the f irst line. 

- the amplitude of two successive lines lies below 50% of the amplitude of the peak and the last line rises again «> 
foot point b line in front of the first line. 

Psak Criteria 

A peak must additionally fulfil the following criteria: 

- minimum peak width = three lines (approx. 0.36 Hz); and 
maximum peak width = 4 Hz. 

For each peak determined in this way. a plurality of specific characteristics will then be determined, which are used 
for the future classification of the said peaks: 
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In the following, the designations used in said table will be explained in detail, 
time: signal time since measurement criteria; 



needle (peak) flfl: 
allTime: 

allFreq: 

fCentr [Hz]: 
fGrav[Hz]: 
reJH[%]: 
absH [norm.]: 



sequence of the peaks/needles found; 

this line does not refer to a peak, but to the whole (preprocessed) time signal in the current FFT 
window; for the regression analysis the respective temporal samples are taken instead of the fre- 
quency lines; 

this line does not refer to a needle, but to the whole spectrum of 0.5 to 10 Hz; for the regression 
analysis, all the frequency lines contained therein are used; 

frequency of the longest line of a peak; 

gravity frequency of all lines of a peak; 

relative height of the peak in comparison with all peaks found; 

absolute height of the peak as normalized after the preprocessing; corresponds to the modula- 
tion degree or approximately to the perfusion; 
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Pslim [factor]: 
ASym [fectr]: 
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slimness factor of the peak, defined as width from foot point to foot point in Hertz; 
measure for the asymmetry of a peak, defined as 



(fo-fa)-ifa'fu) 



w 
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Corel [coeffl: 



75 slope [ratio]: 
Sp02[%]: 

so Perf [%]: 
25 Nadel Score: 



wherein fo = upper foot point, fu = lower foot point, fg = gravity frequency of the peak and w « 
width of the peak; 

correlation coefficient of a needle; i.e. after linear regression analysis of S R , S (R , from the left to 
the right foot point of a peak; describes the slimness of the needle in the x,y representation (Fig. 
4). but not of the peak in the amount spectrum; 

slope of the regression line; corresponds to the ratio; 

Sp02 value determined on the basis of the slope (according to the known empirical relationship, 
cf. equation 2); 

perfusion index determined on the basis of the height of S R » § m* at the point fCentr (e.g. by 
means of the known formula 

Perf = 0.116 • 2S R + 0.626 • 2S IR ; 

number of points according to a classification which will be explained hereinbelGw. 



The peaks, i.e. needles in the x,y representation, which have been identified and characterized in this way in 
accordance with the described embodiment of the method according to the present invention, must now be subjected 
to a selection, since the peaks found will precisely correspond to the fundamental and harmonic waves of the useful 
signal only if the signal is undisturbed. Normally, it must, however, be assumed that in addition to the useful signal com- 
ponents there are arbitrary disturbance components or that arbitrary disturbance components are superimposed on the 
useful signal. 

On the one hand, disturbance components can spectrally be superimposed on the useful signal peak and can 
therefore be an integral component of a needle. On the other hand, disturbance components can be present as sepa- 
rately recognizable needles and peaks, respectively. The first-mentioned case is, in principle, uncritical in the case of 
the method according to the present invention as long as it is assumed that the disturbance within the peak is constant 
to a certain degree. Such a disturbance only results in a broadening of the needle and in a displacement from the origin. 
Unless the correlation coefficient is much smaller than 1 , this type of background disturbance is substantially eliminated 
by the method according to the present invention, i.e. the regression analysis. 

However, in order to eliminate disturbance peaks which appear separately in the distance spectrum, it will be nec- 
essary to subject the peaks, which have been obtained on the basis of the preceding peak identification, to a classifi- 
cation. 

In the following, special relevance criteria that can be used for classifying the peaks will be explained in detail. It is, 
however, apparent that this description only represents a special embodiment; deviating from this embodiment it is also 
possible to use only some of the criteria described, to define the limits for these criteria in a different way, and to award 
points differently. 

First of all, the correlation of the time signals (allTime) must be utilizable, e.g. Corel > 0.4, since otherwise the (per- 
fectly correlated) useful signal component is not sufficiently large; in this case, the signal should not be further utilized 
at all for this time. 

A peak and the needle associated therewith, which are adapted to be used as useful signal, must have all the fol- 
lowing properties: 

the peak fits well into a harmonic frequency series comprising one or several other peaks, and the saturation values 
of the harmonic wave do not differ greatly from one another: 

there are as many harmonic waves as possible; 



the needle is slim, i.e. its correlation coefficient is dose to 1 ; 
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- the frequency of the fundamental wave is in the physiological or specified range (plus tolerance); 
. the saturation value lies in a theoretically useful physiological range (plus tolerance); 

the perfusion lies in the physiological range; and 

- the pulse rate. Le. the peak frequency, lies in the physiologically probable range for the patent monitored, e.g. 
neonate vs. adult 

A peak caused by a disturbance. i.e. a disturber needle, is particularly conspicuous due to the fact that it fails to 
fulfil one or several ones of the above-demanded properties to a special degree. 

The degree of fulfillment of the individual criteria can now be judged in the manner known by a K.0. and a point prin- 
ciple For this purpose, the following K.0. criteria can, for example, be used. i.e. a peak that does not fulfil these criteria 
is not taken into account correlation coefficient > 0.5; pulse rate between 30 and 360 bpm; Sp02 between 0% and 
is 110%; and perfusion index between 0.05 (%) and 25 (%). 

In addition, points can be awarded for the fulfillment of the following criteria: 

- +20 points for each suitable harmonic wave, if the frequency deviation of said harmonic wave is less than 10% (or 
less than 5% for the third and fourth harmonic waves) from the fundamental wave and if the Sp02 deviation thereof 

20 is less than 10% from said fundamental wave; 

+ 10 points for a correlation coefficient which is larger than 0.9; 

- +30 points rf the frequency is in the probable pulse frequency range. 

The peaks and the associated needles can also be evaluated on the basis of tendency variations with regard to a 
reference value that was ascertained in a preceding run. In so doing. e.g. +20 points can be awarded to a respective 
peak for small frequency variations (<1 0% rel). small perfusion variations (<10% rel) or small saturation variations (<5% 
abs) 

in addition, it is also possible to use K.0. criteria for excessively large tendency variations with regard to a reference 
value from the preceding algorithm cycle when the score for a peak and the associated needle, respectively, is compar- 
atively low. Hence, peaks will be rejected, which fulfill the criteria following hereinbelow and which have a lew point 
value. 
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35 - within a short period (15 seconds) after the start of the value output: 

- frequency variation outside of -20% rel ... +30% rel 
. perfusion index variation outside of -20% rel . .. +40% rel 
saturation variation outside of -1 0% abs ... +5% abs; 



in the period following said period (>15 seconds): 



frequency variation outside of -40% rel ... +80% rel 
perfusion index variation outside of -40% rel . .. +80% rel 
45 - saturation variation outside of -15% abs ...+30% abs; 

When a summed needle score has been ascertained in this way for each peak. i.e. for each associated needle, the 
values of the peak whose score is the highest and amounts to at least 60 points are used for outputting in accordance 
with the embodiment described. On the basis of the characteristics for this peak. cf. table 1 hereinbefore, the saturation 
value Sp02 for this peak is then outputted as the saturation value ascertained. In addition, the gravity frequency of this 
peak can be outputted as the pulse rate ascertained and the perfusion index of this peak can be outputted as the per- 
fusion ascertained. ^ , _ . . 

Optionally, it is then possible to filter individual ones or each of the above-mentioned values ascertained. TTie value 
triple can. for example, be inputted together with the last two value triples of the preceding algorithm cycles into a 
median filter so as to eliminate possible outliers, if an incorrect needle should, by mistake, not have been filtered out by 
the above identification and classification method. The median triple ascertained by the median filter then provides the 
reference values which have been described hereinbefore and which can be used for evaluating the tendency. For 
smoothing strongly dynamic or noisy values, an averaging filter. e.g, a 5 value deep box-car averaging filter, can addi- 
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tionally be connected downstream of said median titer. 

If a needle cannot be found by means of the method mentioned hereinbefore, or if none of the needles has a pre- 
determined minimum number of points, the old set of values can be maintained for a predetermined period of time, 
which may be in a range of from 20 to 40 seconds. Subsequently, a special report, INOP report is produced, which 
shows that no utilizable values are produced for the time being. The values can be set to T or "0" in this way. 

A special case can occur if, in spite of a useful signal, no peaks or needles are found because e.g. a very arhythmic 
pulse pattern exists. In this case, a special treatment can be carried out. Although it is impossible to calculate a pulse 
rate due to the arhythmic pulse pattern, an Sp02 value can still be calculated. For use, it is then demanded that the 
correlation coefficient is better than 0.98 in the time domain (allTime) and in the frequency domain (allFreq) and that the 
average value of the saturation values from the alffime and frequency correlation does not deviate more than 5% from 
the reference value. This average value is then outputted as Sp02. 

Fig. 6 shows an overview of a preferred embodiment of a method according to the present invention in the form of 
a flow chart First of all, the raw signals R, IR, which have been obtained from the intensities, are inputted in a process- 
ing unit at 100 as has been described hereinbefore. These raw signals R, IR are subjected to preprocessing at 1 10 so 
as to perform primarily the so-called baseline correction. Subsequently, the values obtained in this way are subjected 
to an FFT, the resultant amplitude spectra for red and infrared being shown at 120. As shown at 130, these amplitude 
spectra are converted into an x,y representation; this corresponds to a complex representation of the amplitude spec- 
tra. A distance spectrum 1 40 is obtained from the amount of the conplex combinatorial values. From this distance spec- 
trum 140, individual peaks, which fulfil relevant criteria, are selected. Subsequently, these peaks are classified e.g. in 
the form of a table 150. On the basis of the characteristics of the individual peaks, said peaks have points awarded 
thereto at 160 so as to identify in this way the peaks that belong to a blood pulse. Following this, a median filtering, an 
average value formation or a historical evaluation 170 are carried out subsequently. Finally, an Sp02 value, a pulse rate 
value and/or a perfusion value are outputted as a result of the above selection, or as a result of the amplitude spectra 
values 120 determined by the above selection, or as a result of the complex combinatorial values 130. rf necessary, out- 
put of an (NOP. instead of an output of the respective values, is effected. 

In the following further developments of and alternatives to the above-described special embodiment of the method 
according to the present invention are deserted. 

It is, for example, possible to determine the above-described peek determination in the red and/or infrared spec- 
trum separately. For this purpose, the ratio following hereinbelow is formed for each sample of the FFT with a (R and a* 
as a Fourier coefficient: ratio(f) = a R (Q/a m (f) . For determining the desired SpQ2, the following frequency components 
are then excluded: all the frequency components for which a Fourier coefficient is smaller than an absolute minimum 
value so as to avoid quantization problems. Furthermore, all frequency components below a relative detector threshold 
are excluded. The detector threshold could be defined as a descending curve, eg. an 1 A shape, 1/f 2 shape, e' x shape, 
in dependence upon the maximum Fourier coefficient a max , e.g. 1 12 a^. Only the frequency values above said thresh- 
old are then used for a ratio Subsequently, a mean ratio is formed by the median, whereupon the standard deviation of 
all ratios is calculated. Above a specific standard deviation, an INOP is outputted instead of an Sp02 value. Following 
this, all ratios which differ from the median by a factor of the standard deviation, e.g. one standard deviation, are elimi- 
nated. The residual ratio points are used for forming the average value. This results in the Sp02 searched for. The pulse 
rate would be the frequency at the maximum value a,^. 

Claims 

1 . A method of determining at least the concentration of a component from the intensity of electromagnetic waves with 
at least two selected wavelengths which are reflected by human tissue or transmitted through human tissue, said 
method comprising the following steps: 

a) converting the intensities of the received electromagnetic signals into at least one first and one second time- 
dependent electric signal (Rf, IRJ; 

b) time-discretely transforming the first and the second electric signal into the frequency domain for determin- 
ing first and second spectral values (S^, S (Rk ) of the first and of the second signal: 

c) forming complex combinatorial values (A^ ration) from said first and second spectra) values: 

d) selecting physiologically relevant combinatorial values by evaluating the complex combinatorial values 
according to given criteria for the physiological relevance thereof: and 

e) calculating the concentration of the component making use of the selected combinatorial values or making 
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use of frequency-wise-associated spectral values. 

2 A method according to claim 1 , wherein, prior to the time-discrete transformation into the frequency domain, said 
first and second time-dependent signals (R|, IRj.) are subjected to preprocessing for removing time-dependent drift 
5 components from said first and second signals. 

3. A method according to daim 1 or 2. wherein the time-discrete transformation is carried out by means of a Fourier 
transformation making use of a suitable time window. 

io 4 A method according to one of the claims 1 to 3 t wherein the selection of the physiologically relevant combinatorial 
values is carried out on the basis of maximum value ranges in the distance spectrum (AO of the complex combina- 
torial values. 

5. A method according to one of the daims 1 to 4 ( wherein the concentration of the component being determined is a 
is gas saturation, particularly the arterial oxygen saturation. 

6. A method according to one of the daims 1 to 5, wherein a pulse rate is additionally determined making use of the 
selective combinatorial values or making use of frequency-associated spectral values. 

7. A method according to one of the daims 1 to 6, wherein a perfusion index is additionally determined making use of 
the selective combinatorial values or making use of frequency-associated spectral values. 

8 A method according to one of the claims 4 to 7, wherein as given criteria for the physiological relevance of the com- 
plex combinatorial values at least two of the following ones are used: width of the maximum value range; frequency 
of the maximum value in the maximum value range; gravity frequency of all combinatorial values in the maximum 
value range; position of the maximum value range with regard to further maximum value ranges in the amount 
spectrum: saturation value obtained for the maximum value in the maximum value range; perfusion index deter- 
mined from said maximum value in said maximum value range; pulse rate determined from the gravity frequency 
in said maximum value range. 

9 A method according to claim 8. wherein steps a) to e) are cyclically repeated, the given criteria comprising in addi- 
tion the following ones: deviation of the frequency of the maximum value of the maximum value range and/or of the 
perfusion index and/or of the saturation value from reference values determined in the course of a preceding cycle. 

10. A method according to one of the daims 1 to 9. wherein steps a) to e) are cydically repeated, the concentration 
" values obtained during several cycles being subjected to a filtering and/or to averaging. 

1 1 A method of determining at least the concentration of a component from the intensity of electromagnetic waves with 
at least two selected wavelengths which are reflected by human tissue or transmitted through human tissue, said 
40 method comprising the following steps: 

converting the intensities of the received electromagnetic signals into at least one first and one second time- 
dependent electric signal (Rj t IRj); 

time-discretely transforming the first and of the second electric signal into the frequency domain for determin- 
ing first and second spectral values (Sr*. S IRk ) of the first and of the second signal; 

selecting physiologically relevant first and second spectral values by evaluating said first and second spectral 
values according to given criteria for the physiological relevance thereof; and 

calculating the concentration of the component making use of the selected first and second spectral values. 

12. A method according to claim 1 1 , wherein the concentration of the component being determined is a gas saturation, 
particularly the arterial oxygen saturation. 

55 

13. A method according to claim 1 1 or 12, wherein, in addition to the concentration of the component, also a pulse rate 
and/or a perfusion index are calculated making use of the selected first and second spectral values. 
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14. An apparatus for determination of at least the concentration of a component from the intensity of electromagnetic 
waves with at least two selected wavelengths which are reflected by human tissue or transmitted through human 
tissue, comprising: 

means for conversion of the intensities of the received electromagnetic signals into at least one first and one 
second time-dependent electric signal 

(Rj, IRj); means for time-discrete transformation of the first and of the second electric signal into the frequency 
domain for determining first and second spectral values (Sp*. S )Rk ) of the first and of the second signal; 

means for formation of complex combinatorial values (A|<, ration from said first and second spectral values; 

means for selection of physiologically relevant combinatorial values by evaluating the complex combinatorial 
values according to given criteria for the physiological relevance thereof; and 

means for calculation of the concentration of the component making use of the selected combinatorial values 
or making use of frequency-wise-associated spectral values. 

15. An apparatus for determination of at least the concentration of a component from the intensity of electromagnetic 
waves with at least two selected wavelengths which are reflected by human tissue or transmitted through human 
tissue, comprising: 

means for conversion of the intensities of the received electromagnetic signals into at least one first and one 
second time-dependent electric signal (Rj, IRj); 

means for time-discrete transformation of the first and of the second electric signal into the frequency domain 
for determining first and second spectral values (Sp^, S (Rk ) of the first and of the second signal; 

means for selection of physiologically relevant first and second spectral values by evaluating said f rst and sec- 
ond spectral values according to given criteria for the physiological relevance thereof; and 

means for calculation of the concentration of the component making use of the selected first and second spec- 
tral values. 
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